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We consider for the first time R-parity violating interactions of the Minimal Standard Supersym- 
metric Model involving neutrinos and quarks ("flavor changing neutral currents", FCNC's) in the 
infall stage of stellar collapse. Our considerations extend to other kinds of flavor changing neu- 
trino reactions as well. We examine non-forward neutrino scattering processes on heavy nuclei and 
free nucleons in the supernova core. This investigation has led to four principal original discover- 
ies/products: (1) first calculation of neutrino flavor changing cross sections for spin 1/2 (e.g. free 
nucleon) and spin nuclear targets; (2) discovery of nuclear mass number squared (A 2 ) coherent 
amplification of neutrino-quark FCNC's; (3) analysis of FCNC-induced alteration of electron cap- 
ture and weak/nuclear equilibrium in the collapsing core; and (4) generalization of the calculated 
cross sections (mentioned in 1) for the case of hot heavy nuclei to be used in collapse/supernova and 
neutrino transport simulations. The scattering processes that we consider allow electron neutrinos 
to change flavor during core collapse, thereby opening holes in the v e sea, which allows electron 
capture to proceed and results in a lower core electron fraction Y e . A lower Y e implies a lower 
homologous core mass, a lower shock energy, and a greater nuclear photo-disintegration burden for 
the shock. In addition, unlike the standard supernova model, the core now could have net muon 
and/or tau lepton numbers. These effects could be significant even for supersymmetric couplings 
below current experimental bounds. 



I. INTRODUCTION 



In this paper we consider the effects of certain Supersymmetry (SUSY) interactions during the core collapse and 
explosion stages of a supernova. These interactions violate lepton number or baryon number and the particular ones 
we are concerned with allow neutrinos and/or antineutrinos to change flavor by scattering off d quarks. We will discuss 
the consequences that such flavor changing interactions might have for current core collapse supernova models. We 
note that our considerations apply to neutrino flavor changing interactions in general, not just those arising from 
SUSY models. 

As the understanding of core collapse supernovae is far from complete, constraints on these SUSY couplings, or 
couplings for flavor changing neutrino interactions in general, do not follow from our considerations. However, we do 
show how such interactions can result in significant alterations in the current model for supernova core collapse and 
explosion. It is conceivable that a better understanding of the supernova phenomena, or the detection of a galactic 
supernova neutrino signal, could allow the effects discussed in this paper to be turned into hard constraints, and/or to 
be used as signatures for R-parity violating SUSY. Likewise the same conclusions would apply to other nonstandard 
flavor changing neutrino interactions. Here, by hard constraints we mean those that can be taken with the same 
reliability and confidence level as constraints derived from terrestrial and accelerator-based experiments. 

R-parity violating interactions have previously been investigated in the context of coherent forward scattering in 
the case of the sunfilp, LJ lJ ll M an d f° r the late stages of supernovae, r-process nucleosynthesis, and the associated 
neutrino signal 0, IE These papers considered quantum mechanical oscillations between neutrino flavor eigenstates. 
The approach in these papers was similar to the Mikheyev-Smirnov-Wolfenstein mechanism|lfj| for matter enhanced 
neutrino mixing due to Standard Model (SM) interactions. 

By contrast, in this paper we look at the interactions in an entirely different context. In particular, we examine 
elastic and inelastic non-forward neutrino flavor changing scattering during the core infall epoch and immediately 
after core bounce. Neutrinos become trapped in the core during infall because of the large number of scatterings 
they have on heavy nuclei. The cross section for elastic coherent scattering of neutrinos on heavy nuclei via the SM 
weak neutral current (NC) is large compared to the cross sections for other types of scattering. Hence, the number 
of NC nuclear scatterings is larger than the number of other types of scatterings. It is this type of scattering which 
is principally responsible for the neutrinos having transport mean free paths which are small compared to the size of 
the core. This regime corresponds to "neutrino trapping." 

The SUSY interactions we consider behave in an analogous fashion to the SM NC neutrino scattering case. In 
particular, the flavor changing SUSY interactions allow for elastic coherent scattering of neutrinos on heavy nuclei. 
We will show that the cross section for this process has much the same form as the cross section for the SM neutral 
current process. However, there are two important differences between these cross sections. The first is that the SUSY 
cross sections are smaller than the SM cross sections because the SUSY interactions are weaker [H E3 than the 
SM interactions. The second difference is that the SUSY interactions allow neutrinos to change flavor. For example, 
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an electron neutrino can change into a muon or tau neutrino by scattering off a nucleus. 

We point out in this paper that even though the interactions are weaker, non forward SUSY scattering events are 
still important because of the large number of scatterings neutrinos undergo in the stellar collapse environment and 
because these interactions allow for flavor changing. This flavor changing is different than the flavor transformation 
(that arises from quantum mechanical neutrino oscillations) which has been studied in the works mentioned above. 
The number of scatterings is large in the post trapping collapse stage because of coherent scattering on heavy nuclei. 
Post bounce, the number of scatterings is large in the neutron star because the density is so high. 

The effects of changing neutrino flavor in the infall epoch of supernova core collapse has been treated in Ref . 01 ■ I* 1 
that work, coherent matter- and neutrino background-enhanced active-active neutrino flavor mixing/transformation 
during the infall epoch was examined. General nonstandard neutrino effects and nonstandard neutrino interactions 
in supernova core collapse are discussed in Ref. [l4| 

Allowing v^s to change flavor can alter the electron capture and thermal physics in the core. In fact, the dynamics 
of core collapse is sensitive to the electron fraction because relativistically degenerate electrons provide most (> 90%) 
of the pressure in the core. By allowing electron neutrinos to change flavor, holes are opened in the Fermi-Dirac sea 
of iVs. These holes allow the electron capture reaction e~ + p <-> n + v e to proceed to the right. The entropy of the 
core, homologous core mass, and the initial shock energy at core bounce all depend on the number of electrons in 
the core during collapse. We will review how decreasing Y e can change these quantities, and what implications such 
changes might have for the supernova model as a whole. We perform a simple analysis to estimate the decrease in Y e 
during infall due to the SUSY flavor changing interactions. 

Another obvious consequence of allowing flavor changing neutrino interactions to proceed during infall is that seas 
of i/„'s and/or z^ r 's with net lepton numbers can be built up. This is in stark contrast to the standard supernova 
model where at core bounce and subsequent to it, there are zero net mu and tau lepton numbers, but a significant 
electron lepton fraction of about Yl = Y e + Y„ e « 0.35. An initial lepton number in mu and tau neutrinos could 
alter all of the subsequent neutrino evolution history in the core and could alter the neutrino signal as well. Lowering 
the trapped electron lepton number fraction post bounce also results in an altered equation of state and changes in 
neutrino transport. We will discuss incorporating these interactions in numerical simulations in order to account for 
feedback on the system and get a more detailed picture of the effects throughout all these regimes. 

The paper is organized as follows. Sec. II will be devoted to SUSY and a derivation of the interactions we are 
considering. In Sec. Ill we discuss the nuclear physics of these interactions and derive relevant cross sections. In Sec. 
IV we discuss the supernova model and give our analysis. In Sec. V we give conclusions. 



II. SUPERSYMMETRY 

Supersvmmetrv[T5l H^ . ITtI ITsj is a proposed symmetry of particle physics. Stated simply, its main feature is that 
each particle has a partner with opposite spin statistics. Stated another way, each fermion particle has a boson 
partner, and each boson has a fermion partner. There are many different models involving SUSY, each containing 
different numbers of unknown parameters. If SUSY is indeed a symmetry of nature, it has not been detected yet. 

Supersymmetry is attractive for theoretical and computational reasons. For example, the underlying symmetry 
group of the SM is the Lorentz group. The underlying symmetry group of SUSY is obtained by adding a generator 
(which represents the super-partner particles) to the generators of the Lorentz group. Thus the underlying symmetry 
group of SUSY is a natural and simple extension of the Lorentz group making SUSY a natural and simple extension 
to the SM. As another example, the three coupling constants of the SM are not equal at low energies, and remain 
distinct when computed at higher energies. In SUSY models, these couplings are still distinct at low energies but 
become equal when calculated at higher energies. In other words, SUSY is a candidate for a Grand Unified Theory. 
Finally, some SM calculations lead to infinities and additional techniques are required to deal with them. However, 
the same calculations in SUSY do not lead to infinities because contributions from particles and their partners cancel 
each other (due to minus signs arising from opposite spin statistics). It is a feature of SUSY that these infinities are 
not present. For this and other reasons, SUSY is taken quite seriously as an extension to the SM. 

The SUSY model we are interested in is the Minimal Standard Supersymmetric Model (MSSM) [3 [HJ with 
additional interactions which break R-parity [2(3, Hi). The MSSM contains the minimum number of fields to describe 
the known 37] SM particles and their superpartners. R-parity is an additional discrete symmetry which requires that 
an interaction must have an even number of SUSY particles. If, for example, R-parity is conserved, then the lightest 
supersymmetric particle (LSP) would not be allowed to decay. This is prevented because the LSP decay would involve 
a vertex with two SM model particles but only one SUSY particle. R-parity violating interactions also violate lepton 
number (L) or baryon number (B). 

The R-parity violating interactions we consider come from adding the following L violating terms to the 
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superpotential 38] : 

A,,/. /.,/.,/-.';. (1) 

Kjk L iQi D k ■ (2) 

Here Li and Qi are lepton and quark SU (2) doublet superfields, E? and £)? are lepton and quark singlet superfields, 
and the A's are coupling constants. (We will later specialize our discussion by considering these couplings to be 
real.) The Roman subscripts are family indices and there is an implicit contraction of SU(2) indices of the doublet 
superfields. Note that in Eq. JQ) A^ is antisymmetric in i and j. 

In general, a superfield consists of scalar, fermion or vector fields and is a function of spacetime coordinates and so- 
called superspace coordinates. The scalar and fermion component fields of a particular superfield represent a particle 
and its super-partner. The superspace coordinates appear as anticommuting Grassman variables. A superfield can 
be expanded in terms of its component fields and these Grassman variables. Likewise, products of superfields such as 
those that appear in Eq.s Q and (J2J can be expanded resulting in products of individual particle and sparticle fields. 
Kinetic and potential terms can be formed from products of superfields and these can also be expanded in terms of 
component fields. A supersymmetric Lagrangian is obtained through a procedure where kinetic, potential and any 
additional interaction terms are expanded and then the superspace coordinates of the result are integrated over 39J. 

Applying this procedure to Eq.s (QJ and J2J gives the following interactions in the Lagrangian pi): 

c = + 44*4 + (4)*^i4 

-(i~j)]+H.c (3) 

l = Ay*i44 + 444 + (4)*^14 
-§i44 - 444 + 

+H.c. (4) 

Here, sparticle scalar fields are denoted by a tilde. For example, <i£,z4,>4 are handed down- type squark, left 
handed neutrino type fermion, and right handed electron type fermion fields, repsectively. Note that where two left 
handed fields are contracted, the first field is charge conjugated. The interactions we are considering in this paper 
come from Eq. (JIJ. In particular, the second and third terms in the first line of Eq. Q each involve a neutrino, 
down type quark and down type squark. These terms and their Hermitian conjugates give flavor changing neutrino 
scattering with d quarks at the tree level. The vertices of these four terms are given in Table [IJ where for illustrative 
purposes, we have taken the couplings to be real. 

We can form tree level diagrams out of the first two vertices in Table [I] Such diagrams are given in Fig.s and 
(0) , where we have chosen exchange of a b squark. Since the indices are unspecified we can change neutrino flavor by 
choosing different indices for the initial and final neutrinos. In the limit of low energy scattering, we can neglect the 
squark momentum in its propagator and then the amplitude for the process in Fig. Q is 

n _ A i31 A £31 J ,.i r i j 

U — 2 d R v L v L dR 

ill.- 

b L 

= -%^^7^i^7M^, (5) 

where the second line has been Fierz transformed. The general low energy effective Lagrangian for our interactions 
can likewise be formed and is given byQ 

A' A' 

r ijk imk £ a i 7m ,j 

2m 2 v L a Ll^ a L 

-%^7^L47m^. (6) 

d L 

The A' coupling constants appearing in these interactions have been constrained with the assumption of 100 GeV as 
a lower bound for all squark masses |40j. 

We are really interested in products of A"s such as those that appear in Eq. ©. In some cases, limits on individual 
A"s are used to constrain products of A"s while in other cases particular products of A"s are constrained. The 
products of A"s for the flavor changing neutrino scatterings we are interested in are constrained to be less than the 
range 10" 2 - 10 5 0,[llL[i3 

The Lagrangian in Eq. © gives neutrino interactions with quarks. In the next section we derive the cross sections 
for neutrino scattering with nuclei and free nucleons. 
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III. NUCLEAR PHYSICS 

In this section we fold the SUSY leptonic and hadronic currents of the Lagrangian in Eq. JBJ into the physics of 
nuclei in a hot medium. In particular, we discuss the coherent elastic scattering cross sections for nuclear target states 
with angular momentum and 1/2 and we present detailed forms for their low momentum transfer limits. We also 
preface this discussion with a brief exposition of the physics of the inelastic SUSY flavor changing (FC) interaction 
channels involving nuclei. 

As pointed out by Bethe, Brown, and co-workers over twenty years ago, the salient feature of the gravitational 
collapse infall epoch of Type II supernova progenitors is that the entropy is low, and remains low throughout collapse. 
In units of Boltzmann's constant k, the entropy per baryon is s ~ 1, some ten times lower than it is in the center of 
the sun. An immediate consequence of this is that nucleons will tend to reside inside large nuclei. Free neutron mass 
fractions will be only of order ~ 10%, while free proton mass fractions will be even smaller in the mildly neutron-rich 
conditions expected in collpase. Therefore, with the rise of density as the collapse proceeds, the mean nuclear mass 
will become larger and larger. 

The neutrino-nucleus cross section in the coherent elastic scattering limit for the Standard Model neutral current is 
proportional to the square of the nuclear mass. Therefore, the trend of increasing nuclear mass eventually causes the 
core's opacity to neutrinos to become large enough for neutrinos to be "trapped." The large opacity causes trapping 
because neutrino mean free paths become less than the physical scale of the collapsing core and likewise mean diffusion 
(random walk) times become longer than the collapse timescale. 

Trapping sets in when the central density of the core is in excess of ptrap > 10 gcm~ 3 . Somewhere between a 
density of 10 13 gem -3 and 10 14 gem -3 , a series of phase transitions will take place where nuclei merge into sheets or 
tubes of material at nuclear density (the so-called "pasta phase"). At higher density, these entities will merge into 
homogeneous nuclear matter. 

The temperature during the infall epoch will likely not change drastically from a value T ss 2 MeV. This is because 
heavy nuclei can store a large amount of energy in excited many-body states, thereby providing the medium with 
what amounts to a large specific heat. We will discuss the general composition and thermodynamic state for the 
stellar core and the prospects for entropy generation in section IV below. 

Nuclear statistical equilibrium (NSE) obtains in the stellar collapse environment, so there may be a fair range of 
nuclear masses represented in the medium. The spins and parities of these target nuclei could span a large range. 
Nevertheless, the two spin cases we consider here in the elastic coherent scattering limit will serve as a guide to the 
behavior of the nuclear cross sections in general and, in particular, encompass the case of free nucleons. Typical 
mean nuclear masses in the density regime where SUSY interactions could play a significant role are A rs 100 to 
200, though the actual NSE distribution of nuclear abundances is sensitive to entropy, neutron excess, the history of 
electron capture reactions and neutrino transport [2~1 125| . 

A complicating issue in estimating neutrino-nucleus cross sections in the inelastic channel is that the target nuclei 
in stellar collapse will not be in their ground states. The large temperatures in the collapsing core, together with the 
large level densities that characterize heavy nuclei, imply that a typical nucleus will be in a state with an excitation 
energy of many tens of MeV above ground, possibly ~ 100 MeV. The mean nuclear excitation energy will be roughly 
the product of the number of nucleons excited above the nuclear Fermi surface (Fermi energy ~ 40 MeV) , and the 
typical energy of each of these excited nucleons. The first of these quantities is ~ aT, where the level density parameter 
is a {A/%) MeV -1 for nuclear mass A. If the nuclear level density is high enough, we can approximate the nucleons 
excited above the Fermi sea to be in plane wave states. In this limit, each excited nucleon has an energy of order the 
temperature, so that the mean excitation energy is (E) f=a aT 2 . 

Neutrino-nucleon interactions in collapse will nearly always involve a nucleon inside a nucleus. The amplitudes 
for neutrino-nucleon interactions, whether mediated by SM or non-standard processes, will therefore likely involve 
nontrivial medium effects on account of the dense environment where the nucleons and nuclei are found. However, 
we can follow the usual procedure for Standard Model nuclear weak interactions in collapse and assume that nucleons 
in nuclei will have neutrino interaction properties broadly similar to those of bare nucleons, but restrict ourselves to 
processes that are kinematically allowed in the medium. To this end, we derive neutrino-nuclear elastic scattering 
cross sections for nuclei and bare nucleons. 

We derive the cross sections in terms of a general low energy effective Lagrangian for neutrino-quark interactions. 
We define what the general operators in this Lagrangian are for i) the Standard Model and ii) our SUSY model. We 
then give the cross sections in terms of the general operators and comment on these specific cases. 

Throughout, operators are implicitly evaluated at x = 0. For example, q^^q stands for g(0)7 M g(0). We take particle 
states to be relativistically normalized, 

{p'\p)=2E{2n) 3 S^{p-p'). (7) 
We will use k, k' for neutrino momenta and p,p' for the momenta of nuclei. Here, unprimed and primed variables are 
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for incoming and outgoing particles, respectively. We also recall here the invariant Mandelstam variables s, t and u 
(defined by N(p) + v(k) -► N(p') + v(k')) 



s = (p + fc) 2 = (p' + fc') 2 
i = (p- p ') 2 = (fc-fc') 2 



with 



s + 1 + u = 2_. '• 



(8) 
(9) 
(10) 



In the rest frame of the nucleus s, t, and u are 



s = M Z + 2ME V 



t = 



-AM El sin 2 (6/2) 
M + 2E V sin 2 (6»/2) 



9 4Af£ 2 sm 2 (<9/2) 
u = M 2 -2E U M 1 



Af + 2£„sin 2 (6 , /2)' 

where 6* is the scattering angle, M is the nuclear rest mass and the neutrino mass has been neglected. 



(11) 
(12) 

(13) 



A. Effective Lagrangian 

We can write a general expression for the low energy effective Lagrangian for neutrino-quark interactions as 



£ e G r = [i?V(l - 7sK] EfcVWl- + K% l5 )q], 



(14) 



where the sum runs over all quarks. For the Standard Model (SM) neutral current, the Lagrangian 26] has the usual 
current-current form 

Gf 



C cS = -J= [P e7 "(l - 7>e] 

x2 ^ g 7M (T 3 - sin 2 0„Q) g 



(15) 



where the sum runs over left and right handed helicity states of all quarks. For this SM case, the operators in Eq. 
(fHjl are, 



K 



(q) _ GF (rr ; ., 2 , 



— ^(T 3 -2sin 2 ^Q)<5, 



.(9) 



— T 3 <5 
V2 



(16) 
(17) 



In Eq.s (|16|l and (|17fl . T 3 and Q are the isospin and charge of quark q, respectively. For example, we have T 3 = 
l/2(-l/2) and Q = 2/3(-l/3) for g = u(<2). 

For the R-breaking SUSY Lagrangian (in particular, that given in Eq. 10 but where the couplings are complex) 
the operators as defined in Eq. (|14f> are 



8 ^ 



Kik(^'jik)* Kki(^'jki)* 




(18) 
(19) 
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In summary, the basic current-current structure of the low energy effective Lagrangians for the Standard Model 
flavor-conserving and SUSY flavor-changing neutrino interactions are identical. This affords a simple extension of the 
usual nuclear systematics of the Standard Model interactions to putative SUSY-inspired flavor changing reactions. 
It should be kept in mind that other non-SUSY flavor changing neutrino-quark interactions may have a completely 
different form with, consequently nuclear response characteristics quite different from those presented in the following 
subsection. However, our derivation of the coherent elastic cross sections in the subsequent subsections will apply to 
any flavor changing interactions that can be described by the Lagrangian in Eq. 1141) . 



We first discuss general (i.e., not necessarily coherent elastic) neutrino scattering from nuclear states. We will 
exploit here the similarity between the basic SUSY neutrino interaction operators and the Standard Model ones, 
as was discussed above. Again, we warn the reader that other, non-SUSY inspired flavor changing neutrino-quark 
interactions may have quite different nuclear operators. 

We describe the ith excited state of a nucleus with Z protons and N neutrons with a general many-body state 
ket \ 1 $i(p)) which can be thought of as a superposition of Dirac spinors for individual nucleons. Here the nuclear 
wave function ^ describes the full many-body nuclear state which, in the low energy nonrelativistic limit would have 
good quantum numbers corresponding to excitation energy Ei, angular momentum/parity J? , and isospin Ti with 
z-projection Tf = (Z — N)/2. The kincmatical bulk center of mass momentum of the state is labeled by p. 

Of course, given the temperature and density conditions expected in stellar collapse as outlined above, and given 
the relatively modest neutrino energies expected to be encountered, we are not likely, at least initially at neutrino 
trapping, to require the full relativistic description of nucleons in the nuclei. Therefore, we will begin our discussion 
of the nuclear cross sections by considering general (e.g., inelastic) neutrino scattering on nonrelativistic nuclei and 
nucleons. 

The small components of the full nuclear many-body wave functions/spinors can be taken to be linear superpositions 
of Slater determinants of nonrelativistic two-component single particle wave functions. These single particle states we 
envision to be computed in the usual independent particle model regime (e.g. shell model or RPA), as Hartree-Fock 
solutions with a specified two-body nucleon-nucleon potential. Likewise, the coefficients in the linear superposition of 
Slater determinants are presumed to come from an overall diagonalization procedure involving the residual interaction 
between nucleon quasi-particles in the nucleus. 

The formal nuclear currents for general neutrino-nucleus scattering in the fully relativistic regime for either the 
Standard Model or SUSY-inspired neutrino interactions are: 



Here Vt/ is the many-body field operator describing a final nuclear state obtainable by the scattering channel being 
considered. This could differ from the initial nuclear state by a difference in any of the quantum numbers mentioned 
above. 

In particular, a change in excitation energy between the initial and final nuclear states constitutes an inelastic 
neutrino interaction. This could, in turn, be classified as either of two types. It could be an endothcrmic transition, 
where Ef > Ei, and the scattering neutrino gives up energy to the nucleus. It could be an exothermic transition, 
where Ef < Ei and the scattering neutrino picks up energy from the (excited) nucleus. The latter process is possible 
because, as discussed above, the target nuclei of interest in stellar collapse are expected to be in very highly excited 
states. 

If we take the nonrelativistic limit for the nucleon spinors in Eq.s (|2U[1 and (|21[1 . and contract with the leptonic 
current with an assumed negligible momentum transfer (the allowed approximation), then we obtain the usual Fermi 
and Gamow- Teller forms for the nuclear matrix elements, the absolute squares of which are 



B. General Nuclear Matrix Elements 



V = */(i/)7 M *i(p) 
A" = %(p')7 M 75*i(p). 



(20) 
(21) 




(23) 



(22) 



Here "0i and ij)f are the initial and final two component nuclear many body wave functions which correspond to 
the nonrelativistic reductions of \& j and 'J / , respectively. In the above definitions we have averaged over initial and 
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summed over final states for the particular transition i — > /. The sum on n is a sum over nucleons in the nucleus and 
r is an operator which depends on the process mediating the transition. The total matrix element squared between 
initial and final nuclear states i and / can be written 

|A%| 2 = |D F | 2 |Mf/ + \D GT \ 2 \M^f, (24) 

where Dp and Z?gt are coupling constants, each of which depends on the particular scattering mode, as does the 
operator r. 

For example, in the Standard Model charged current weak interaction we have Dp = Cy, the vector coupling 
constant, while -Dgt = Ca, the axial vector coupling constant. In this case r is the isospin raising or lowering 
operator for individual nucleons, so that ^ n r(n) — T , the overall isospin raising/lowering operator. The Fermi 
matrix element between members of an isospin muliplet in the Standard Model charged current case is simply 

\M? f f = KV/IT^)! 2 (25) 
= T(T+l)-T 3 (T 3 ±l) . (26) 

The selection rules implied by this are AJ = 0, and no parity or isospin change. By contrast, the Gamow- Teller 
operator for the Standard Model case, J2 n T ( n )&> ^ s a spatial vector and an isovector, so that the selection rules are 
A J = 0, 1 (but no — > 0), no parity change and AT = 0, 1 (no — > 0). We discuss the Standard Model neutral 
current analogs of these interactions below in our presentation of the coherent elastic scattering case. 

The SUSY flavor changing neutrino interactions are a different mode, with different values for the couplings and 
different meaning for the operator r. We shall denote the couplings in the SUSY flavor changing mode as Dp — Cy USY 
and Dgt = Cqt SY (which can be expressed in terms of the unkown couplings, A"s ). For the <i-quark flavor changing 
neutrino scattering channels discussed above, the operator r simply counts d-quarks. In this case then, r(n) = 3/2+t z , 
where t z is the z-component of isospin for individual nucleons. Note that this operator yields r(n) = 2 for neutrons 
and 1 for protons. The Fermi matrix element between members of an isospin multiplet is 

\M? f f = 2N + Z, (27) 

with selection rules the same as those for the Standard Model Fermi operator. Likewise, the selection rules for the 
Gamow- Teller operator in this case are the same as in the Standard Model case. 

In both the Standard Model and SUSY cases the form of the operators determines selection rules and the general 
distribution of strength with excitation energy. Note that the Fermi operator commutes with the nuclear Hamiltonian 
if we neglect the Coulomb potential as small compared to the strong interaction. On the other hand, the Gamow- 
Teller operator does not commute with the nuclear Hamiltonian. This is because of the spin dependence of the nuclear 
Hamiltonian. As a result, the "Gamow- Teller strength" associated with a given initial nuclear state in general will 
be spread throughout the other states of the nucleus that meet the selection rules. In practice, it is known that in 
the Standard Model Gamow- Teller case the strength may be significantly collected in a narrow resonance. We would 
expect this for the SUSY case as well, but the different dependence of the SUSY operator on isopsin probably implies 
that the resonance excitation energy centroids, widths, and overall strength will be different. 

For larger momentum transfer the operators and the corresponding matrix elements will depend on the momentum. 
For SM weak interactions, the momentum transfer dependence of the one-body operators and, hence, all nuclear 
matrix elements is classified by a comparison of inverse momentum transfer to nuclear size R. The "allowed" regime 
corresponds to values of this parameter \ftR < 1, while "forbidden" weak interactions have larger values. For typical 
pre-bounce collapse conditions, neutrino energies are E v ~ 50 MeV, while nuclear masses are large A > 100 (with 
nuclear radii given by R ~ 1.2A 1 / 3 fm). The high Fermi levels of either the e~ or v e sea (and sometimes the internal 
nuclear weak selection rules) dictate relatively low momentum transfer. With modest momentum transfers, the 
allowed approximation may be reasonable for both the electron capture and ^-capture channels of the SM charged 
current weak interactions. (See Ref. [24| for an in-depth discussion of these points.) 

Flavor changing scatterings involving initial state f e 's likely will be somewhat different than in the SM neutral 
current case for nuclei. Momentum transfer in this channel could be larger on account of little or no blocking of 
final state and/or v r phase space. This difference is accentuated for inelastic and endothermic neutrino scattering 
channels, where the nuclear final state is at a higher excitation energy than the initial state, and the neutrino final 
state has a consequently lower energy than the initial state. For flavor-preserving v e interactions, the final state v e is 
apt to be blocked by the degenerate v e sea unless its energy is larger than the v e Fermi level. This is less likely to be 
the case for flavor changing interactions because there are essentially no or v T species in the infalling supernova 
core, absent flavor changing neutral currents (FCNC's). 

As the collapse proceeds and densities and neutrino energies rise, momentum dependence in the operators and 
matrix elements may become important. This would tend to skew the effects we discuss in the Analysis section to 
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the highest energy i/ e 's if the SUSY couplings are small, but at the highest densities or larger coupling it will make 
little difference. 



C. Spin-0 Coherent Elastic Nuclear Scattering 

We now consider elastic scattering where the initial and final nuclear states are the same. We here designate the 
ket describing a spin-0 nuclear state with total momentum p as \(Z, N),p). This state is some linear combination of 
the states for individual quarks and nucleons which make up the nucleus, as discussed above. We need the matrix 
element of the hadronic current in Eq. I|14|l with respect to such a spin-0 nuclear state. The matrix elements for the 
vector and axial vector pieces expressed in terms of form factors are: 

((Z,N)p'\q^q\(Z,N)p) = f q (t)(P+P'Y, (28) 
((Z,N)p'\q^ l5 q\(Z,N)p) = 0, (29) 

where t is the momentum transfer squared defined in Eq. i|12|) ■ The matrix element of the axial vector current is zero 
because the initial and final states have the same parity and thus the left side of Eq. (|29|l is odd under parity, but 
the vectors p^ and p' M cannot be combined on the right hand side into an object that is odd under parity. 
The form factor f q (0) is fixed. The conserved charge jUl associated with the vector current is 

Q q = J d 3 xqry°q. (30) 

For example, this charge for the u quark is Q u = J d?xvJu. The spin-0 nuclear states given above are eigenstates of 
the charge operators defined by Eq. (J3UJ), and in particular, these charges count number of quarks. Therefore, setting 
fi = in Eq. (|28|l and using Eq. 0, we have f q {0) = N q , where N q is the number of g-quarks in the (Z, 7V)-nucleus: 

/ tt (0) = 2Z + N (31) 
f d (0) = Z + 2N. (32) 

The precise momentum dependence of the form factor in Eq. Ij28(l is not known of course, but we may be able to get 
insights into this by considering the SM case in the context of typical stellar collapse conditions, as discussed above. 
In accordance with the arguments made there about the role of forbidden transitions in collapse, in what follows we 
will take f(t) ~ /(0). As we will see below, this is well justified in the case of nucleons, for which the form factors 
are precisely known and the approximation of treating them as constant introduces an error no larger than 3.0%. 

For the general low energy effective Lagrangian given in Eq. (|14() in the channel Vi + (Z, N) Vj + (Z, N) and/or 
V{ + (Z, N) — > Vj + (Z, N) the differential cross section is 



datj 1 

dt 7T 



4«/«i(*)+«Vii/u(*)| S 



ts 



(s - M 2 ) 2 



(33) 



Integrating over squared momentum transfer t, the total cross section is 



.^Kijfdio) + /4y/«(o)r — ^ (34) 

2 '4i j fd(°)+ K vijfu(0)\ 2 El (35) 



7T 



Note that for the case of the SM neutral current (with couplings given by Ea.sll6land ll7fl this result agrees with the 
cross sections given in Ref. 27] . For the R-parity violating SUSY model, where we have only a d-quark interaction, 
the term Ky^fu is not present in the preceding equations. In particular, the cross sections for this case, with operators 
given by Eq. (|18[1. are 



l'K Vij \ 2 f d (QfEi. (36) 

The results in Eq.s l|33|) - (|35|) would apply as well to any general neutrino flavor changing interactions that can be 
described by the Lagrangian in Eq. (|14fl . 
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D. Spin-| Nuclei 

For a target nucleus with initial state angular momentum J = 1/2, the matrix elements expressed in terms of form 
factors are: 



((Z,N)p's'\q^q\(Z,N)ps) = u»V 



u s \p') 



((Z,N)p's'\q 7 ^ 5 q\(Z,N)ps) = u s \p') 



Hi{t)^ + ^i(jp-p') v a^ 



2M 



{Hi{t)-H 3 {t)yt 



u s (p), 



2M 



(p + p'T 

l5U S (p). 



(37) 
(38) 
(39) 



The form factors in this case are H\{t), H^it), G\(t) and Ga(i). The states are labeled by their momentum p, spin 
s, and u s (p) is a Dirac spinor satisfying the Dirac equation (7 ■ p — M)u s (p) = 0. In going from Eq. I|37|) to Eq. I|38|) 
we used the Gordon identity. Eq. (|37|l is the standard form, but Eq. 138|) is simpler for cross section computations. 
Hence, our results are given in terms of the combination 



H 1 = Hi — H2 



(40) 



The form factors are implicitly labeled by the flavor of the quark in the current. Note that from the Dirac equation 
and conservation of momentum we have 



(41) 



In computing the cross section, we will neglect neutrino masses m v . As a result of this approximation, Gi will be 
irrelevent. 

Current conservation gives 



#i(0) 



2Z + N for q = u 
Z + 2N for q = d 



(42) 



If we replace the vector current q'y^q by the electromagnetic current J^ m = ^ii'y^u — ^d^^d, the form factors Hi 2 are 
replaced by the well known electromagnetic form factors Pi,2- In that case, we have more information. In particular, 
we know in this case that the magnetic moment in units of e/2M is fi — -Fi(O) — F 2 (0). 

Moreover, for nucleons we have two more pieces of information. They form an iso-doublet and the electric and 
magnetic factors, 



Ge = Fl ~ 4M^ 
Gm = Fi — F 2 , 



F 2 



are known empirically to have common g 2 -dependence, 



G EM {q 2 ) = G EM (0)g(<l 2 ), 

m = 1 



[l-q^/Mlf 



(43) 
(44) 

(45) 
(46) 



where M 2 » 0.71 GeV 2 . 

So, in the case of nucleons, we know precisely the form factors for both the dj^d and the vr/^u currents. We define 
H^ n) by 



(n\ d^d In) = (p\ wy^u \p) — u s (p 1 ) 
(p\ dj^d \p) = (n| u7 M u |n) = u s (p 1 ) 



r"'''m-/< . g a" (*) 



2M 



i(p-p') v a^ 



H { i p \t)r + ^^ i (p- P 'u^ 



u s (p), 



(47) 
(48) 
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and, as before take Hi = Hi — H 2 . We have used isospin symmetry, which gives two relations among the four matrix 
elements. Using the normalization conditions in Eq. 142II . the electromagnetic form factors and isospin, we solve for 
the individual quark currents. We find 



H 



I (*) 



(2/i p + Mn )5(t) a (3.7)5(t) 
0*p + 2fx n )g(t) « (-l.O)p(i) 
(2 - 2/x p - M „) 5 (t) « (-1.7) 5 (t) 
ff a W (t) = (1 - M P - 2/*n)ff(*) « (2.0) fl (t), 



if. 



where 



1 



(49) 
(50) 
(51) 
(52) 

(53) 



-4E 2 

y.max 



t/AM 2 ' 

It is worth pointing out that H^ p ' n \t) varies little over the range of interest, that is, from t min 
—0.01 GeV 2 to t max = 0. Indeed, (g(Q) — <?(tmin))/<?(0) — 3.0%. To the extent that nucleon form factors are generic, 
the approximation of treating nuclear form factors as constant seems reasonable. 

Now we turn our attention to the axial current. As mentioned after Eq. I|41|). G 2 will not be needed. Therefore, 
we only address G\. For protons we have [28| 

for q = u 



( 0.78 ±0.03 
Gi(0) = <^ -0.48 ± 0.03 for q = d 
I -0.14 ±0.03 fovq = s 



(54) 



while for neutrons the same expressions apply but with u 
factor for the difference 



d exchanged. It should be noted that the t « form 



(p's' \ U7 p 7 5 u - d^^ z d \ps) = u s (p')gA 



7 + — (p-p r 



'(p)- 



(55) 



is very well known from neutron beta decay: gA = 1.2670 ± 0.0035 

The differential cross section for spin-^ nuclei for a general Lagrangian which has only neutrino interactions with 
with d quarks is, 



\ K Vij\ 2 



2fff (1 - t/AM 2 ) + 4H 2 Hi [(a - M 2 ) 2 + st] + H\ [2(s - M 2 ) 2 + 2st + t 2 ] 



daij 1 
~dT ~ 2vr(s - M 2 ) 2 ' 

+ \kmj\ 2 G 2 i [2(s - M 2 ) 2 + 2st + t 2 - 4M 2 t] ± 2Re(^ ljK %)HiGi [2t(s - M 2 ) + t 2 ] 



(56) 



The ± takes the upper sign for neutrino scattering and the lower sign for antineutrino scattering. The corresponding 
total cross section, approximating H{t) sa H(Q) is then 



2\2 



(s - M 2 ) 
6tts 3 



( jip H 2 s [( s - M 2 ? + l2M 2 s] + 6s 2 HiH 2 + H 2 [As 2 - 2sM 2 + M 4 ] 



1 



+ K 



Aij 



2 G\ [As 2 + AsM 2 + M 4 ] T 2Re( K ^4* J )fi r 1 G 1 (s - M 2 ){M 2 + 2s) 



(57) 




QHiH 2 + mi) +9|«^ ij |2 



G\ 



(58) 



(The accuracy of these cross sections in the free nucleon channel may be improved by utilizing the well determined 
form factors of Eq.s 49 - 52.) Note that the cross term {Ky i jK A * i j) is of order E% and hence was neglected in the last 
line. However, it should be kept in mind that this cross term may be important for generating asymmetries. 

For Ky A given by Eq.s l|18fl and (|19l) . the above equations give the differential and total cross section for the R-parity 
violating SUSY model. The results for the SM are obtained by exchanging Ky — > (k v + k v ) and k a — > (k a + k a ) 
everywhere in the above equations and using the appropriate values for the form factors (that is, taking account of 
whether they came from a u quark or d quark current and including cross terms between them). 
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IV. ANALYSIS: FCNC EFFECTS IN CORE COLLAPSE 
A. The Supernova Model 

Gravitational collapse and the transformation of gravitational binding energy into (mostly) neutrinos and a small 
amount of outgoing kinetic energy and radiation is thought to be the process which powers supernovae of Type II, lb, 
and Ic. The energy resident in the degenerate seas of neutrinos in these objects is huge, constituting some 10% of the 
rest mass of the compact object produced in the collapse. In dramatic contrast, the energy of the explosion (optical 
plus outgoing kinetic energy) is only 1% of this. We clearly see that even small changes in the energy or flavor content 
of the degenerate electron and v e reservoir may affect the physics of collapse and explosion. 

The gravitational collapse phenomenon can be exquisitely sensitive to lepton number violating physics. This is 
because stellar collapse generates prodigious degenerate seas of electrons and electron neutrinos in weak equilibrium 
in a well-ordered, low entropy state. These electron lepton number fermion seas dominate the pressure and, hence, the 
dynamics during the collapse process. Flavor changing interactions can convert some of this degenerate electron lepton 
number into seas of neutrinos with mu and tau lepton numbers, thereby altering the Fermi levels of the electrons and 
v e 's and, consequently, changing the pressure and, to a lesser extent, the entropy in the collapsing star. 

Better theoretical and observational insight into core collapse supernovae could allow them to become the ultimate 
laboratories for studying neutrino physics beyond the Standard Model. Conceivably, some day we could obtain 
constraints on or even discoveries of new physics, otherwise unobtainable in conventional accelerator experiments and 
terrestrial laboratories. We are, unfortunately, not currently at this level of understanding. However, this is not to 
say that new neutrino physics could not significantly alter our picture for how stars collapse and explode. 

The current paradigm for the core collapse supernova explosion mechanism involves gravitational collapse halted 
by a "bounce" at nuclear or super-nuclear density. The bounce of an inner core is accompanied by the generation of 
an initially energetic shock wave at its boundary. This shock has its energy degraded by nuclear photo-disintegration 
of heavy nuclei as it traverses the outer core (the material laying between the inner core surface and the edge of the 
initial iron core). As a result, the shock is weakened and becomes a standing accretion shock, incapable of exploding 
the star. It is thought, however, that this shock subsequently is revived and strengthened by neutrino heating, perhaps 
aided by convective and hydrodynamic processes |29tl30l ]. 

The progenitor star of such an event would have mass M > 10 M Q . This star would evolve over some millions 
of years through a succession of nuclear burning stages, eventually producing an "onion skin" structure of layers of 
fossil ashes of each core burning stage along with active burning shells at their boundaries. At the center would be a 
core composed of iron peak material in Nuclear Statistical Equilibrium (NSE). The central density of this core will 
be p ~ 10 10 gcm -3 , or pio ~ 1, and the central temperature will be T ~ 1 MeV. The electrons in these conditions 
are relativistically degenerate and they supply nearly all of the support pressure. The electron Fermi energy is 

1/3 

/i e w 11.1 MeV(/9ioi / e ) , where Y e is the electron fraction or net number of electrons per baryon. Comparing /i e to 
the temperature T suggests that the entropy is low. A detailed accounting of all degrees of freedom yields an entropy 
per baryon s w 1 in units of Boltzmann's constant 23] . 

The silicon burning shell will add iron peak material to the core, until it exceeds the Chandrasekhar mass (1.2 Mq 
to 1.6 Mq in this case), whereupon the core will go dynamically unstable and collapse at an appreciable fraction of 
the free fall rate. As the collapse proceeds, the baryons principally will reside in heavy nuclei, as dictated by the low 
entropy and as described in detail above in section III. Forcing the baryons to be confined to heavy nuclei implies 
that the baryonic pressure contribution during collapse is nearly negligible. 

During this collapse process, electron capture on protons which reside in these big nuclei 

e- + A(Z,N) ^ A{Z -l,N + l) + v e (59) 

and on the very few free protons present will lower the core electron fraction Y e and produce v e 's. Electron capture 
on nuclei tends to increase the entropy because it leaves daughter nuclei in states that lie above the mean thermal 
excitation energy by ~ 3 MeV, a typical spin-orbit nuclear shell splitting for representative nuclei early in the collapse. 
The entropy increase via this process is modest and is nowhere near large enough to melt the nuclei. Only completely 
melting the nuclei would allow the baryonic pressure to become significant. In turn, this would require an entropy 
increase of at least three units of Boltzman's constant per barvon [l3llT^ |. 

The collapse is initially homologous. That is, the infall velocity of a fluid clement is proportional to its radius 
from the center. The collapse time scale is of order the free fall time, in particular, a few seconds. As the collapse 
proceeds and the pressure is reduced (because electron capture reduces Y e ), only a smaller, inner core can continue 
to collapse homologously. This inner core region of mass 0.6 — 0.8M© remains homologous, while the outer region of 
mass 0.7 — 0.9M© collapses supersonically. 

The edge of the homologous core is roughly the sonic point in the infall velocity field. In a sense, the inner core 
is causally self connected via pressure (sound) waves. Therefore, this portion of the core will bounce as a unit 



12 



when the central density becomes large enough for the nuclei to merge and the pressure becomes dominated by the 
nonrelativistic baryonic component. This will occur, typically, at 3 to 5 times nuclear matter saturation density. 
Another consequence of the causal structure of the inner core is that its mass is essentially that of an instantaneous 
Chandrasekhar mass, 

M hc « 5.8F e 2 M . (60) 

This is a good approximation because electron degeneracy pressure dominates. 

As outlined in section III above, the growth of large nuclei during core infall eventually causes the ^ e 's produced by 
electron capture to become trapped in the core. This occurs when their random walk, or diffusion times become longer 
than the collapse time scale. The principal arbiter of trapping initially is coherent elastic neutral current neutrino 
scattering on the heavy nuclei. This process has a cross section proportional to the square of the nuclear mass and 
the square of the neutrino energy. As a consequence, at first only the higher energy v e 's are trapped. Lower energy 
neutrinos can still escape and take lepton number and entropy out of the core. 

Neutral current coherent neutrino scattering also is conservative. It is only when non-conservative neutrino-electron 
scattering and neutrino- neutrino scattering opacity sources become significant, and when hot nucleus de-excitation 
into neutrino pairs operates that the v e distribution evolves toward a thermal, equilibrated Fermi-Dirac form. In this 
case the core approaches weak, or beta equilibrium. At this point neutrinos are truly trapped and the evolution of the 
core from this point on is nearly adiabatic. Note that neutrino pair (y a , v a with a =e,/i,r) production via electron 
bremsstrahlung and plasmon decay processes during collapse are suppressed on account of the low entropy and the 
extreme electron degeneracy and high Fermi level fi e . As a result, during infall there may be relatively fewer neutrino 
types other than v e in the medium. 

In weak equilibrium the electron capture and reverse v e capture reaction rates are equal and larger than the collapse 
rate. In essence, the build up of a degenerate v e sea eventually blocks further electron capture. Any process that 
opened holes in the v e distribution would inevitably lead to further electron capture and a lower Y e and, hence, a 
smaller homologous core mass at bounce, with the deleterious consequences outlined below. 

In weak equilibrium the chemical potentials of leptons and nucleons are related through 

Me - Hu e ~ A + Sm np , (61) 

where (j, = /2n — /i p is the difference of the kinetic chemical potentials (not including rest mass) of neutrons and protons 
and Sm np ~ 1.293 MeV is the neutron-proton rest mass difference. In NSE, the nucleon chemical potentials are the 
same inside and outside nuclei. Inside nuclei /2 can be interpreted as the difference of the neutron and proton nuclear 
Fermi levels. This rises as electron capture proceeds and the medium and the nuclei become more neutron-rich. The 

1 /3 

electron neutrino Fermi level (chemical potential) is roughly fx Ue w 11.1 MeV(2pxoY Ue ) , where the net number of 
u e 's over P e 's per baryon is Y Ve . 

Simulations of hydrostatic stellar evolution coupled with weak interaction rates suggest that the electron fraction 
at the onset of collapse is Y™ lt w 0.42. Simple estimates of self consistent electron capture and nuclear equation of 
state issues yield a typical electron fraction at neutrino trapping of about Y" e tra P « 0.35. Subsequently, lepton capture 
reactions redistribute this electron lepton number so that when weak equilibrium obtains shortly before core bounce 
Y e « 0.3 and Y Ve w 0.05, to give a trapped total electron lepton number fraction equal to Y* rap . 

At a time of about 10—100 milliseconds after the neutrinos are trapped and attain beta equilibrium, the inner region 
of the core reaches nuclear density and halts its collapse. The supersonically-infalling outer core material bounces 
off the inner core and a shock wave forms at their boundary. The shock's initial energy can be approximated as the 
kinetic energy of the infalling outer core material. This kinetic energy, in turn, is approximately the gravitational 
potential energy of the inner core. In terms of Y e in the core at bounce, the initial shock energy depends on the 
electron fraction roughly like [24| 

E^L ~ (n) 10/3 . (62) 

The initial shock energy is ~ 10 51 ergs in the standard collapse model, a figure tantalizingly close to the energy of 
the shock observed to emerge through the progenitor star's envelope. However, entropy considerations again dictate 
nuclear physics realities that result in the vitiation of the initial shock strength. 

All of the infall kinetic energy of the inner homologous core is converted to heat at core bounce. Furthermore, at 
core bounce all of the energy stored in the excited states of heavy nuclei (see section III) is returned to the general 
medium of homogeneous nuclear matter created from the merged nuclei. The upshot is a high temperature for the 
shocked core, T « 20 MeV to 70 MeV|3lj. 

The entropy jump across the shock front is about a factor of ten. In these hot, disordered conditions behind the 
shock, NSE favors disintegration of nuclei, so-called nuclear photo-disintegration. However, since each nucleon is 
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bound in a nucleus by on average 8 MeV, the shock loses energy as it propagates through the low entropy material 
in the outer core. In fact, the shock will lose ~ 10 51 ergs for each 0.1 M Q of outer core material traversed. The shock 
propagates through the outer core and the material beyond. Not surprisingly, simulations show that the shock stalls, 
or more accurately becomes a standing accretion shock, a few hundred km from the center of the core. This all 
happens within a few milliseconds after core bounce. 

In this condition the shock cannot effect an ejection of the outer stellar material to space. That is, it cannot 
explode the star. However, there is still a huge resevoir of neutrino energy in the core. Some of this energy could be 
transported and deposited behind the shock, reviving it and leading to an eventual explosion. 

The higher entropy and high temperatures of the post-shock, post-bounce core (the hot proto-neutron star) now 
allow rapid production of neutrino pairs of all flavors. The neutrino and antineutrino seas for muon and tau neutrinos 
are approximately the same with no net lepton number in cither one, but there remains an excess of z/ e 's over D e 's on 
account of the trapped electron lepton number in the star. 

There is rough equi-partition of energy among the six active neutrino species (y e ,v e ,v ll , : D lli ,v T ,v T ') in the hot proto- 
neutron star. These neutrino species diffuse out of the core on a random walk time scale Tj iff ~ seconds, so that 
the energy luminosity in each of the six neutrino species is a few times 10 52 ergss _1 during the first half second or 
so and falls to ~ 10 51 ergss _1 some 10 s after bounce. These represent huge neutrino fluxes. In fact, the goal of 
a successful supernova explosion theory is to explain how 1% of this energy can be pumped into the shock, either 
through direct heating of the material behind the shock, perhaps aided by convection, or through alleviation of some 
of the photo-disintegration burden on the shock. 

In summary, the initial iron core with a radius of order the size of the earth collapses to a proto-neutron star with a 
radius of ~ 40 km in about a second, promptly liberating ~ 10 52 ergs of energy, 90% of which resides in neutrino seas 
of all flavors, and 10% of which, some 10 51 ergs, comprises infall kinetic energy which becomes the initial shock energy. 
The shocked core, now a hot proto-neutron star, undergoes subsequent quasi-static contraction on a neutrino diffusion 
(random walk) time scale, T&s ~ 3 s. On this time scale the hot proto-neutron star will contract to a radius of about 
10 km and release 10 53 ergs in gravitational binding energy. More than 99% this energy will appear as neutrinos of all 
kinds. The basic points of this scenario have been confirmed by the neutrino-induced events detected from SN 1987A. 

Flavor changing neutral current (FCNC) neutrino processes can affect the above scenario for gravitational collapse 
and shock generation and propagation significantly. If electron neutrinos change flavor then holes are opened in their 
Fermi Dirac sea allowing electron capture to proceed. Therefore, the number of electrons will be reduced and, hence, 
less pressure will be provided by electrons. The mu and tau neutrinos created from electron capture plus FCNC 
processes will contribute some extra pressure but, as a result of their much smaller Fermi levels (chemical potentials), 
this will not compensate for the reduction of electron pressure. 

A lower electron fraction and consequently lowered pressure on infall results in a double whammy. Not only is the 
initial shock energy reduced as Eq. (|62|) shows, but the homologous core mass is reduced, leaving a larger outer core 
and a consequently greater nuclear photo-dissociation burden for the shock to overcome. In the popular current model 
for core collapse supernovae, this would greatly lessen the chance of getting a viable shock and an observationally 
acceptable explosion. 

And there is another consequence of the operation of flavor changing neutrino processes during core infall. The 
FCNC processes could allow the post-bounce core to possess significant net mu and tau lepton numbers. This could, 
in turn, alter the hot proto-neutron star de-leptonization history through feedback on the nuclear equation of state 
and neutrino transport. This would also alter the flavor content of the emergent neutrino fluxes and therefore possibly 
change the supernova neutrino signal. 

B. Evaluation of Cross Sections 

In this section we use the cross sections derived previously to estimate the number of flavor changing (FC) scatterings 
which neutrinos might undergo at different densities during core collapse. These estimates are not meant as accurate 
predictions for the actual number of scatterings because we do not account for feedback in the system. Instead, 
these estimates are meant to illustrate the possible consequences of FCNC interactions and to suggest where a proper 
treatment through numerical simulation is warranted. For large enough SUSY couplings, we could have enough 
flavor changing events during collapse and in the post bounce neutron star to result in a core evolution and shock 
generation/propagation scenario significantly different from that in the standard supernova model with only SM 
neutrino interactions. To illustrate these results, we limit our discussion to the case of spin zero nuclei with Z = N. 
This case is arguably overly simplistic, but it serves to illustrate how and at what level new flavor changing neutrino 
processes can be important. 

The cross sections we need are given in Eq. (|35H for the SM neutral current and in Eq. (|36[) for the SUSY model. 



14 



The SM cross section for the case of Z = N can be written as 



(63) 



where A = Z + N and 



_ 4G 2 F m 2 e h 2 



1.76 x 10" 44 cm 2 . 



(64) 



Note that G 2 F has been absorbed into ctq. The SUSY model cross section can likewise be written as 



713 ~ 4 G 2 



[2N + ZYI — 
1 m P 



(65) 



where we have defined 



V2 



8 



"4 



(66) 



Note the similarity between between G l FC and the Fermi constant which is defined as 



Gf 

71 



8m 2 / 



(67) 



Both are of the same form: a dimensionless coupling constant squared divided by a mass squared. It would be 
especially desirable to define a flavor changing SUSY coupling Gfc = |A'|/8m|. This could be factored out of Eq.s 
(|18fl and (|19|) . and likewise Eq. (|66|) . Then, any differences in the A"s might be described by some underlying 
symmetry for the model. However, as the coupling constants are unknown we will use the notation of Eq. I|66(l to 
stress that we are dealing with flavor changing neutrino interactions that behave like the SM interactions. 

To evaluate the SUSY cross sections we need values for the constants G % FC defined in Eq. 1)66(1 . Values for G FC 
are obtained by substituting values for individual A"s, or specific products of A"s, with indices corresponding to the 
particular neutrino flavors being considered. 

The A' coupling constants are constrained from many different experiments [TH Il2l EH] . As can be seen from Eq. 
a particular A' with one set of indices serves as the coupling constant for several vertices. Thus, several different 
processes can involve the same couplings and more than one experiment can be used to constrain a particular coupling. 
Hence, in the literature, claimed constraints on the couplings can vary. Some constraints from Ref. [TlT ] are: 



|Ai lfc | < 2 x 1CT 2 <-> 2.9 x lCT 1 ^ 



R 

1 ~fc 



|Ai fel |<4x 10- 2 ^7.1 xlO- 1 ^ 
\\' 21k \ <6x 10- 



1.5 x ictM 



-13k 
R 



lA^Kl.SxlO-X 



\\' 3lk \ < 1.2 x lO- 1 ^ 

The notation |A'| < n x d k is shorthand for |A'| < n x (mjfc/lOOGeV 



consistent with the constraints from references P, 0, y, 



7. 



(68) 
(69) 
(70) 
(71) 
(72) 

We point out that these constraints are 
The A"s above can be combined to give 



G FC for 



Vt 



flavor changing processes. If we ignore phases and assume that no drastic cancellations occur in Eq. (|66|l . then for 
the values of A"s given above, the upper limit on G l p C can range from 10- 6 GeV~ 2 to 10^ 8 GeV" 2 El- 

As discussed above, the dominant neutrino scattering opacity source after trapping is coherent nuclear scattering. 
The number of scatterings is then given approximately by the random walk expression 



R 
at 



(73) 
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where R is the radius of the core and At is the distance a neutrino travels before scattering off target nucleus T, i.e., 
the mean free path between such scattering events. The neutrino mean free path after neutrino trapping and before 
the formation of the nuclear "pasta phase" is given by 

(A T ) _1 ps n T a T , (74) 

where nx is the number density of target nuclei T, and or is the cross section for coherent scattering off such nuclei. 
The number of flavor changing (FC) scatterings 7V s ^ tt that neutrinos undergo after trapping can be estimated from 
the number of neutral current (NC) scatterings. Substituting Eq. (|74fl into Eq. (|75|) and taking the ratio N^ tt /N^ tt 
gives 

N FC ( f^4°\ N NC /yrN 

Jv scatt ~ ^NC J JV scatti \ 10 ) 

where ct^ c and a^ c are the FC and NC scattering cross sections on a nuclues of mass number A. As alluded to above, 
to illustrate the effect of FC scatterings we only need to consider neutrino scattering off nuclei with Z — N and so we 
use the cross sections given by Eq.s and <|HH|l . Since (27V + Z) 2 ~ A 2 , Eq. lf75^l with these cross sections becomes 

N ¥C ( G FC \ N NC (76) 

V G Z F Sill' 9 W J 

Taking the appropriate value of the Weinberg angle from experiment gives sin 2 6 W = 0.23 which, together with the 
Fermi constant G-p = 1.166 x 10 _5 GeV -2 , allows Eq. l (75|l to become 

Nl c att « 1 x 10 21 [G FC GeV 2 ] 4 i\££t- (77) 

With this expression we can calculate the number of FC scatterings as a function of the number of NC scattering 
for various values of G l p C . We note that similar estimates for the number of scatterings for general flavor changing 
neutrino interactions could be obtained in this same way once they have been cast into the nuclear physics framework 
discussed in section III and their cross sections are known. 

To get the number of NC scatterings we use Eq.s (|73|) and l|74|l . We take the cross section in Eq. (|63[) with a mean 
neutrino energy E v figured from the density in the core. Likewise ha, the number density of target nuclei of mass A, 
can be computed from the core density. The mean free path is |3^| 

(A^ c )- 1 ^4xl0- 5 ^ 3 cm- 1 , (78) 

where p±% — p/10 12 gcm~ 3 . 

We estimate the number of FC scattering from Eq. I|77|) and calculate the mean free path for neutrinos to undergo 
an FC scattering using Eq. I|73|l . For these calcuations we approximate the radius of the core as 50 km over the whole 
range of densities we consider. Keep in mind, however, that the core radius will be some 45 km immediately after core 
bounce, and will shrink over several neutrino diffusion times T^ iff to « 10 km. The rate of FC scatterings is obtained 
here by dividing the speed of light by the mean free path. 

We do these calculations for a range of densities covering the post neutrino-trapping infall epoch of collapse and 
the post bounce hot proto-neutron star. We have used the cross sections in Eq.s l)63f) and l|65|l over the entire range of 
these densities. This is not strictly valid at larger densitites where, as discussed in section III, nuclei do not exist and 
where neutrino opacities may be dominated by scattering on collective modes in nuclear matter or different phases 
of quark matter which can exist at high density. However, coherent neutrino scattering effects with quarks can still 
occur during these phases j34[. Hence, our procedure may give a crude estimate for the number of scatterings during 
these phases, and the level at which FCNC couplings may be large enough to be of significance. 

For our calculations, we will use G FC values of lO^GeV" 2 , lO^GeV -2 , 10~ 8 GeV~ 2 and lO^GcV" 2 . The case 
where we take Gfc = 10 _9 GeV -2 covers the situation where at least one of the A"s is an order of magnitude smaller 
than the smallest value suggested by current experimental bounds. We give our results in Tables HTllVI These results 
are meant only as a guide to the order of magnitude of the number of FC neutrino scatterings and, hence, to the level 
of FC effects in the core. 

The estimated total number of scatterings for a given density is not meant to be the total number for all the 
densities at that order of magnitude. Rather, they are a guide to the total number of scatterings during each of 
several time intervals (which typically are a few milliseconds) as the central density of the core passes through each 
regime. Additionally, results given for higher densities were not computed with values corresponding to a core in 
which FC processes had been operative. In other words, we have neglected feedback here and any semblance of a 
self consistent collapse history when FC effects are appreciable. In our conclusions we discuss more about possible 
additional feedback the FC processes could have on the model and stress again that the numbers in these tables are 
merely a guide. 
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C. Analysis 



The numbers provided in the tables above give insight into the range of possible physical effects of the FC interactions 
on collapse and shock generation/propagation physics and the rough level of FC coupling where these effects can be 
expected to be important. They suggest that for some plausible as yet unconstrained ranges of these couplings, v^s 
can change flavor. Comparing the rates for electron capture (essentially given by the inverse of the weak equilibration 
time scale) and possible FC processes, we see that it is likely that the former process will be fast compared to the 
latter one. In this limiting case, every time a hole is opened in the u e sea by an FC interaction, it is immediately 
filled by a v e created by electron capture. This approximation is used and discussed in Ref. fl3l ]. 

Another limiting case is one in which the FC processes are fast and equilibrium among flavors is achieved and all 
the Fermi levels of the neutrino seas are the same, though not necessarily at the maximum level. This case likely 
would have to correspond to larger values of Gfc- 

FC interactions will tend to change v e to tT and thereby build up the muon and tau neutrino seas, but not vice 
versa. This is because the v e sea is already in weak equilibrium and possesses a rather high Fermi level. Because of 
this significant v e Fermi level, ZAj )T — * v e tends to be blocked. Of course, it is possible that a hole in the v e sea could 
get filled by a u e created in a v^ T FC interaction before it gets filled by a v e created via electron capture. However, as 

a e A cap is larger than <x^ c , holes are more likely to be filled by electron capture. The main point is that after trapping, 
in the limiting case where electron capture rates can be regarded as fast comparted to FC rates, the level of the v & 
sea remains the same as it would be absent FC interactions. 

We can quantify the effect of the FC interactions for this limiting case by using the approximation that the level of 
the v e sea remains the same. The lepton fraction (total lepton number of all kinds per baryon) in the core is given by 

Y L = Y e + Y v ,+Y Vlt +Y Vr , (79) 

where the net number fraction relative to baryons is defined as discussed in subsection IIV Al by. 

rif — tit 

Yf = (80) 

n b 

We can evaluate this at neutrino trapping and again at some time past trapping. At trapping, there is no excess of 
j^.t's over ^ ljT 's because these species have all been produced as neutrino- antineutrino pairs via thermal emission 
processes. Therefore, by Eq. ffity we have 

Y*™ p = F£ ap = 0. (81) 

At later times, if our FC interactions exist and are taking place, then Y v Vt increases while Y Vc remains the same and 
Y e decreases. Therefore, the post-trapping v e fraction is fixed in this limit, 

F pt = y£ a p. (82) 

As discussed above above, after trapping the total number of relativistic leptons is fixed if we neglect neutrino radiation 
from the surface of the collapsing core. Therefore, Yl is constant in this limit and we have the relation, 

y e trap + r£ ap = y pt + rj e rap + y% + y* . (83) 

In other words, the change in electron fraction is given by 

AY e = - (V pt + Yg) . (84) 



The effect of FC interactions essentially is to convert AY e to Y v „ T . 

We do not give specific values for AY e because we are not equipped to account for feedback. We can however point 
out limiting cases as illustrated in Figure |21 The extreme limiting case on infall is 

Yv, = Y Vt = l£ app . (85) 

The muon and tau neutrino seas are built up to the maximum level before bounce. This limit is actually consistent with 
the current constraints on the A"s. For the result of Eq. 1)85(1 to obtain, we require a total number of electron neutrinos 
equal to 2 x A*' app to change flavor. Given that the rate of FC interactions for the value of Gfc = 1 X 10~ 7 GeV -2 
is of order 10 scattering per millisecond for p\2 = 10 we have more than enough room for this situation. For tighter 
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constraints on the A"s which give G-pc of order 10 8 GeV 2 , 1Z ~ 10 scatt/ms at pi2 = 100 with enough time prior 
to bounce to build up the making the extreme limiting case possible. 

For this case, we would expect AY e w —2Y Ve w —0.1. This is a large effect, on the same order as the change 
in electron fraction from the introduction of new electron capture physics |35|. It would reduce the initial shock 
energy by a factor of ~ 3 and increase the nuclear photo-disintegration burden on the shock by a substantial amount. 
This would, in turn, imply that the shock would stall at a smaller radius, likely decreasing the efficacy of neutrino 
re-heating. However, although these effects argue in the standard picture against a viable explosion, caution is called 
for and detailed shock propogation physics must be employed. (Again see Ref. |35) on this issue where it was pointed 
out that consistent electron capture physics on infall and during shock reheating is necessary to gauge the viablity 
of the shock.) Ultimately, we do not understand where explosions come from, so that constraints on SUSY or other 
non-SM FCNC interactions cannot be drawn with any degree of confidence yet. Clearly, to accurately assess the 
effects of FCNC's on collapse/explosion, our cross sections should be incorporated into detailed simulations with full 
hydrodynamics, nuclear equation of state, and Boltzman neutrino transport. 

V. CONCLUSIONS 

We have shown that new flavor changing neutral current (FCNC) neutrino interactions could have significant effects 
in the physics of core collapse supernovae. In particular, we have outlined how these interactions can both reduce 
the initial shock energy and increase the amount of material that the shock must traverse and photo-disintegrate. 
Both of these effects go in the direction of causing the shock to weaken and, therefore, to stall and become a standing 
accretion shock at a smaller radius. This likely would greatly decrease the prospects for obtaining a viable shock and 
explosion in the current paradigm for the core collapse supernova explosion mechanism. 

Another consequence of a significant level of FCNC neutrino interactions during collapse would be the creation 
of net mu and tau lepton numbers in the post-bounce hot proto-neutron star core. This differs from the standard 
supernova collapse scenario, where these lepton numbers would be zero. This could affect neutrino transport and the 
equation of state in the de-leptonizing core. It could also alter the neutrino radiation emergent from the core and, 
thereby change the neutrino signal expected from a core collapse event. 

However, as emphasized at the end of the last section, only detailed collapse simulations with coupled nuclear 
physics, hydrodynamics, and neutrino transport with our FCNC cross sections can adequately ascertain what hap- 
pens in the supernova. We have demonstrated that there are as yet unconstrained ranges of R-parity violating SUSY 
FCNC parameters that could yield alterations in weak interaction physics and therefore the FCNC's are worthy of 
incorportation into these simulations. Likewise, other kinds of flavor changing (non-SUSY) neutrino-nucleus inter- 
actions would produce alterations in lepton physics during collapse identical in a qualitative sense to what we have 
considered here. 

There are also effects which we have not analyzed. For example, anti-neutrino flavor changing due to these in- 
teractions would have similar cross sections, and in particular A 2 -like dependence. Another effect is the potential 
for entropy increase due to flavor changing inelastic scattering. That is, entropy increase due to high energy i/ e 's 
scattering into low energy i^'s (because these low energy states are not blocked) leaving nuclei in excited states and 
thus heating the system. These effects treated properly in a simulation might lead to further interesting results. 

Additional effects we have not analyzed include FC interactions involving the charged leptons. Note that the last 
two terms in the Lagrangian in Eq. and their Hermitian conjugates can mediate — * fjT . During the early stages 
of the post-neutrino trapping infall regime this process will not be operative because the electron Fermi energy fi e 
will be less than the muon rest mass (m^ « 106 MeV). However, near core bounce and post core bounce, where the 
density is of order nuclear saturation density, it is possible to have /i e > m M . In this case we would expect e~ — > 
and for large enough FC coupling equilibrium between these species would result. This would augment the trends 
discussed above which stem from neutrino FC reactions. Namely, with large scale conversion in the channel e~ — > fi~ 
we would expect a further lowering of Y e and increase in net muon lepton number. 

We have outlined the nuclear physics of FCNC reactions. We have found that, similarly to the SM NC case, there 
is a coherent elastic scattering neutrino FC channel with a cross section that scales as (27V + Z) 2 . Given the low 
entropy, and consequently the large neutron-rich nuclei expected in the infall epoch of stellar collapse, this factor can 
be sizable. 

As we have seen, the square of the cross section comes into the calculation for number of scatterings. Since nuclei 
would be even more neutron rich because of increased e~ capture if FCNCs were taken into account, the factor of 2 
in 2N + Z could be important. We neglected this factor of 2 (factor of 4 in the cross sections) and related feedback 
in our estimates for numbers of scatterings. If we included this factor, the numbers of scatterings could increase by 
up to a factor of 16. 

Our considerations for the results of neutrino flavor changing apply to any general FC interactions which give large 
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enough numbers of FC scatterings. Large cross sections would likely arise from neutrino quark interactions that 
allowed for coherent scattering amplifications. For a general interaction with neutrinos on electrons, for example, the 
number of scatterings would likely not be large . Note that our R-parity violating SUSY model has neutrino-electron 
flavor changing interactions. In particular, the second and third terms in the first line of Eq. @ (which originates 
from the superpotential term in Eq. ^) lead to neutrino-electron flavor changing. The products of A's that appear 
here again result in ranges of values for a low energy effective coefficient similar to the ranges we used for Gfc- 

In the case where FCNCs are the result of our R-parity violating SUSY model, even couplings an order of magnitude 
and possibly more below experimental bounds could have significant effects on the infall epoch of stellar collapse. This 
naturally begs the question of whether we can use the considerations presented here to provide constraints on these 
couplings or even, conceivably, to suggest signatures of non-SM neutrino interactions in, for example, the supernova 
neutrino signal. 

We argue that the current lack of understanding of the origin of explosions in core collapse supernovae precludes 
constraints and signatures which can be treated with the same confidence as those obtained from current accelerator- 
based and terrestrial laboratory experiments. However, there is intense interest in core collapse supernovae and new 
insights are being gained from large scale simulations with state of the art numerical hydrodynamics, nuclear physics, 
and neutrino transport |36j . Ultimately, a better understanding of these events may result and this, in turn, could 
allow the considerations discussed in this paper to be turned into legitimate constraints and/or signatures. The 
prospects for this would be greatly increased by the detection of a neutrino burst signal from a collapse event in our 
galaxy. 

On the other hand, the experimental discovery of supersymmetry would necessitate a new assesment of the stellar 
collapse problem. The discovery of R-parity violating SUSY in particular would be a strong motivation for FCNCs. 
In this case the stellar collapse effects of FCNCs outlined in this paper would become potentially vital (depending on 
the strength of coupling constants) in understanding the core collapse supernova phenomenon. 

We would like to point out a cautionary historical tale for the supernova modeling community. In the 1970's the 
discovery of SM neutral currents in the laboratory (and their theoretical prediction) completely altered the model 
for core collapse at the time. Likewise, the existence of new FCNC interactions could significantly alter the current 
model for collapse and explosion of massive stars. 
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FIG. 1: Flavor changing neutrino scattering with a d quark. This tree level diagram is obtained by joining the first two vertices 
in Table Time advances from left to right and the labels are for fermion fields. 




FIG. 2: Flavor changing anti-neutrino scattering with a d quark. This tree level diagram is obtained by joining the first two 
vertices in Table Time advances from left to right and the labels are for fermion fields. 
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Standard Moderate Max 



FIG. 3: Possible scenarios for the effect of FC interactions on the levels of the various fermion seas. The bin with verticle bars 
represntes Y e , that with horizontal bars represents Y Vl! , and those with diagonal bars represent Y v T . "Standard" refers to the 
standard collapse model in weak equilibrium. "Max" refers to the maximal effects of FCNCs as discussed in the text, while 
"moderate" refers to an intermediate case. 
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TABLE I: Vertices of interactions. The first and third terms appear in Eq. (QJ while the second and fourth are Hermitian 
conjugates of the first and third. The indices are assigned as e, d = 1; fj,, s = 2; r, b = 3. The couplings here are taken to be 
real. 
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P12 

10 12 g/cm 3 


N NC 
1 * scatt 


N FC 

Jv scatt 


A^ C = fl/(JV™ t t ) 1/2 
m 


TZ = c/\ 
scatt /ms 


1 


4 x 10 4 


4 x 10 1 


8 x 10 3 


40 


10 


9 x 10 7 


9 x 10 4 


2 x 10 2 


2 x 10 3 


100 


2 x 10 11 


2 x 10 8 


4 


8 x 10 4 


1000 


4 x 10 14 


4 x 10 11 


8 x 10~ 2 


4 x 10 6 



TABLE II: Values for FC interactions for Gfc = 1 x 10 6 GeV 2 Mean free paths are given in meters, rates are given as 
scatterings per millisecond. 



Pl2 

10 12 g/cm 3 


n nc 

Jv scatt 


N FC 

Jv scatt 


A^ c = R/(N™ t )^ 
m 


TZ = c/\ 
scatt / ms 


1 


4 x 10 4 


4 x 10~ 3 


8 x 10 5 


.4 


10 


9 x 10 7 


9 


2 x 10 4 


20 


100 


2 x 10 11 


2 x 10 4 


4 x 10 2 


8 x 10 2 


1000 


4 x 10 14 


4 x 10 7 


8 


4 x 10 4 



TABLE III: Values for FC interactions for Gfc = 1 x 10 7 GeV 2 Mean free paths are given in meters, rates are given as 
scatterings per millisecond. 



Pl2 

10 12 g/cm 3 


n nc 

Jv scatt 


N FC 

Jv scatt 


A^ C = fl/(JV£tt) 1/a 
m 


TZ = c/\ 
scatt / ms 


1 


4 x 10 4 


4 x 10~ 7 


8 x 10 7 


4 x 10~ 3 


10 


9 x 10 7 


9 x 10~ 4 


2 x 10 6 


2 x 10" 1 


100 


2 x 10 11 


2 


4 x 10 4 


8 


1000 


4 x 10 14 


4 x 10 3 


8 x 10 2 


4 x 10 2 



TABLE IV: Values for FC interactions for Gfc = 1 x 10 8 GeV 2 Mean free paths are given in meters, rates are given as 
scatterings per millisecond. 



P12 

10 12 g/cm 3 


/V NC 

Jv scatt 


N FC 

Jv scatt 


Ai c = R/iN™^ 2 
m 


TZ = c/X 
scatt /ms 


1 


4 x 10 4 


4 x 10~ n 


8 x 10 9 


4 x 10~ 5 


10 


9 x 10 7 


9 x 10 -8 


2 x 10 8 


2 x 10 -3 


100 


2 x 10 11 


2 x 10~ 2 


4 x 10 6 


8 x 10~ 2 


1000 


4 x 10 14 


4 x 10 _1 


8 x 10 4 


4 



TABLE V: Values for FC interactions for Gfc = 1 x 10 9 GeV 2 Mean free paths are given in meters, rates are given as 
scatterings per millisecond. 
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A concise discussion of these and other terms which can be added to the superpotential, and implications of adding all 
these terms, is given in Ref. |22^|. 

For more details on superfields and constructing a Lagrangian from a superfield action we refer the reader to Ref.s 
[Til ITU IT^I . For a quicker review we recommend Ref. Il9) . 
This has been done in Ref. |2l| and for a more recent review see Ref. [TH|. 
Note that this charge is not the electric charge. 

Constraints in these papers are given in terms of a parameter e, not on products of A''s directly. Constraints on products of 
A''s can be solved for in terms of e assuming a squark mass of 100 GeV. Note that these constraints are on specific products 
of A''s, not on individual A''s. Note also there are notational inconsistencies in Ref. E| in Table I and the definition of the 
parameter e. 

[43] Note that these values for the coefficients are also consistent with constraints for general FC neutrino interactions given 
in Ref. EH. 



